Dynamic critical exponent of two-, three-, and four-dimensional XY models with 
relaxational and resistively shunted junction dynamics 
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The dynamic critical exponent z is determined numerically for the d-dimensional XY model 
(d — 2, 3, and 4) subject to relaxational dynamics and resistively shunted junction dynamics. We 
investigate both the equilibrium fluctuation and the relaxation behavior from nonequilibrium to- 
wards equilibrium, using the flnite-size scaling method. The resulting values of z are shown to 
depend on the boundary conditions used, the periodic boundary condition, and fluctuating twist 
boundary condition (FTBC), which implies that the different treatments of the boundary in some 
cases give rise to different critical dynamics. It is also found that the equilibrium scaling and the 
approach to equilibrium scaling for the the same boundary condition do not always give the same 
value of 2. The FTBC in conjunction with the finite-size scaling of the linear resistance for both 
type of dynamics yields values of z consistent with expectations for superfluids and superconductors: 
z = 2, 3/2, and 2 for d = 2, 3, and 4, respectively. 



PACS numbers: 74.40. -Hk, 05.70.Jk, 75.40. Gb, 75.40.Mg 



I. INTRODUCTION 



Superconducting films, Josephson junction arrays, and 
superfluid "^He are systems where topological defects play 
an important role close to the phase transition. This 
is particularly striking in two dimensions (2D) where a 
phase transition of the Kosterlitz-Thouless (KT) nature 
is driven by the unbinding of thermaJJi^ created topolog- 
ical defects, vortex-antivortex pairs.tra In 3D such topo- 
logical defects take the form of vortex loops and it has 
been argued that the physics clpse to the transition can 
be associated with these loops.H The common feature in 
these systems is that they can be characterized by a com- 
plex order parameter. The XY model can be viewed as a 
discretized version of such systems where only the phase 
of the complex order parameter plays a significant role. 
This model is believed to catch the essential features of 
the topological defects present in ''He as well as in super- 
conductors in the limit when the magnetic penetration 
length is much larger than the correlation length; high-Tc 
superconductors fall into this category.a All the systems 
which can be described by the XY model belong to the 
same universality class for the thermodynamic critical 
properties of the phase transition. 

In the present paper we have the connection between 
the XY model and superfluid and superconducting sys- 
tems in mind. However, the XY model per se can equally 
well be viewed as a simple model of a ferromagnet where 
the phase angle corresponds to the direction of a 2D spin 
vector associated with each lattice site. 

Our interest in the present paper is the dynamical 
properties associated with topological defects, which may 
of course depend on the explicit choice of the dynamics 
imposed on the model. We here investigate two types of 
dynamics: One is a simple relaxational dynamics (RD) 
and the other is the resistively shunted junction dynam- 
ics (RSJD). We calculate the dynamic critical exponent z 



using various scaling relations both associated with equi- 
librium and with the approach to equilibrium when start- 
ing from a nonequilibrium configuration. Our main con- 
clusion is that the dynamic critical exponents associated 
with the topological defects are the same for these two 
types of dynamics, RD and RSJD. However, this conclu- 
sion does depend on the precise treatment of the bound- 
ary. We demonstrate that various values of z can be 
obtained by changing the treatment of the boundary, as 
well as by changing from scaling in equilibrium to scaling 
for the approach to equilibrium. 

This paper is organized as follows: In Sec. |^ we briefly 
introduce the XY model and explain how the dynamic 
equations are defined in RSJD and RD taking boundary 
conditions into account. Section 
scaling relations used to obtain 

simulations are given in Sec. rtVl for spatial dimensions 
d = 2, 3, and 4, whereas Sec.^ contains discussions of 
the results. Finally Sec. ^ gives a short summary of the 
main conclusions. 



II. XY MODEL AND DYNAMICS 

A. XY model 

The c?-dimensional XY Hamiltonian on a hypercubic 
lattice of the size il = L'^ is defined by 
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where the summation is over nearest neighboring pairs, 
6r is the phase of the complex order parameter at posi- 
tion r, and J is the coupling strength. The XY Hamilto- 
nian is appropriate not only to describe the overdamped 
Josephson junctions arrays without charging energy, but 
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can also be viewed as a discretized form of the Ginzburg- 
Landau (GL) free energy 

FGL[V(r)] =Jdr (^al^(r)p + ^|^(r)p + i|V^(r)p) , 

(2) 

where the amphtude fluctuations of the complex order 
parameter ip{r) are neglected: ?A(r) = ■f/^oe*®'^'"-' with tpo 
fixed to a constant. When mapping the GL free en- 
ergy functional onto the XY Hamiltonian the coupling 
strength J is found to be proportional to jV'ol^- 

The thermodynamic properties of the XY model have 
been intensely studied for many years and it is well known 
that the important length scale in the critical region, 
the correlation length ^, diverges at the critical temper- 
ature Tc- In 3D and 4D the divergence is of the stan- 
dard form of the continuous second-order transition, i.e., 
f(r) ~ ir-Tcl"", whereas in 2D In^(r) - {T-T^)-'^/'^ 
as Tc is approached from above and ^ = oo in the whole 
low-temperature phase where quasi-long-range order ex- 
its in the absence of true long-range order.EJ From the 
point of view of the finite-size scaling, this feature of the 
2D KT transition turns the finite system size L into the 
relevant length scale in the low-temperature phase. 



B. Boundary condition 

Experiments on superconductors and ^He are usu- 
ally done on samples with open boundaries. From this 
perspective it is preferable to use boundary conditions, 
which reflects this experimental situation also in the sim- 
ulations. However, simulations of the XY model can 
usually only be well converged on relatively small lattice 
sizes, and since the surface to volume ratio is inversely 
proportional to the linear system size L, the open bound- 
ary gives rise to large surface effects, which decay very 
slowly as the system size is increased. The standard way 
of reducing these unwanted surface effects is to impose 
the periodic boundary condition (PBC): Or+Lji — ^r, 
where (i denotes the basis vectors of the lattice, e.g., 
jl — x,y, z in 3D. One drawback of this boundary condi- 
tion is that it restricts the twist from r to r + Lfi, defined 
as the sum of the phase differences along a direct path 
connecting the two positions, to an integer multiple of 
2tt. On the other hand, this twist from one boundary to 
the opposite for an open system can have any value. It 
is thus preferable to relax the PBC so as to allow for a 
continuous twist by changing the boundary condition to 
a more generalized form: ^rHi^pi = + L/S.^^, which has 
been used in various contextsJjis In particular the bound- 
ary condition where the twist variable is not fixed to 
a constant but allowed to fluctuate has been termed the 
fluctuating twist boundary condition (FTBC), which was 
originally introduced for static Monte Carlo (MC) sim- 
ulationa3 and then extended to Langevin-type dynamics 



at finite temperatures. la Since the FTBC allows for any 
value of the twist, it is closer to the open boundary con- 
dition for a real system. Of course one does not expect 
the treatment of the boundary to affect the results in 
the thermodynamic limit. However, as we will show and 
discuss here, the dynamics at criticality can depend on 
the boundary condition, in as far as the dynamic critical 
exponent can be defined in terms of the finite-size scal- 
ing. It is worth mentioning that a similar observation, 
i.e., that an important exponent may depend on bound- 
ary conditions, has been made recently in tke study of 
the stiffness exponent of vortex-glass models. □ 

C. Dynamic models 

Next we introduce two simple dynamic models widely 
used to describe behaviors of superfiuids, superconduct- 
ing films, regular Josephson junction arrays, and also 
bulk high Tc superconductors close to the transition tem- 
perature. 

1. Resistively shunted junction dynamics 

A d-dimensional hypercubic array of size = L"^ {L = 
linear size) of superconducting grains weakly coupled by 
resistively shunted Josephson junctions is effectively de- 
scribed by the XY Hamiltonian (|^) when it comes to 
the static properties. On the other hand, dynamic equa- 
tions of motion for the corresponding overdamped RSJ 
model are generated from local conservation of the cur- 
rent on each grain. The total current Irr' between neigh- 
boring grains (r, r') is the sum of the supercurrent, the 
normal resistive current, and the thermal noise current: 
^rr' = -^rr' + ^rr' + ^rr' ■ The supcrcurrent is given by 
the Josephson current-phase relation 7^^., — /c sin((/)rr')i 
where Ic = 2eJ/h is the critical current for a single junc- 
tion. The normal resistive current 7^^., = V^-r' /Ro, where 
the voltage difference T4r' is related to the phase differ- 
ence by Vrr' = {Ti/2e)(j)rr' and i?o is the shunt resistance. 
Finally the thermal noise currents 1*^., in the shunts sat- 
isfy (/*^,) = and 

Uq 

(3) 

where (• • ■) is the thermal average, and 6(t) and ^rr' are 
the Dirac and Kronecker deltas, respectively. From local 
current conservation we obtain 

J2lrr+n=ir, (4) 

n 

where the h summation is over 2c? nearest neighbors of 
site r on a hypercubic lattice in d dimensions (n = ±/i), 
e.g., n — ±x, ±y, ±i in 3D, and I^^^ is an external cur- 
rent source at r (in the present work, we only consider 
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the case without external driving: I^^^ = 0). Introducing 
the lattice Green's function Urr', which is the inverse of 
the discrete Laplacian, the RSJD equations of motion in 
the absence of external currents can be written in dimen- 
sionless form as 



dt 



E 



sin(6'r 



hn) + Cr 



(5) 



where U^r' = U^r' — U^r ,0 and from here on we normalize 
the time, the current, the distance, the energy, and the 
temperature in units of h/2eRoIc, Ic, the lattice spac- 
ing a, J, and J/ks, respectively. The on-site noise term 
Cr(i) = - Y.r'^^^'Y.n^l'v'+hi^) Spatially correlated, 
which is a consequence of the local current conservation, 
and satisfies (Cr(t)Cr' (0)) = 2TUrr'S{t). The RSJD equa- 
tions (|^) can be rewritten in a Langevin-type forir£3 



'dt 



' S9r'{t) 



Cr 



(6) 



with the XY Hamiltonian H in Eq. (|l|) [compare with 
Eq. ( |l^ ) for relaxational dynamics]. 

We now introduce the FTBC for the RSJD. The global 
twist LA jj in the fi direction across the whole system (see 
Sec. p Bp is introduced through the local transformation 
9 J. Or + r ■ A, still keeping Or = Or+Lfi as the periodic 
part of the phases. The Hamiltonian in terms of these 
variables is 



H[Or, A] = 



E 

rp. 



COS fr — 



'r+p. 



(7) 



where the /t summation is over d nearest neighbors of 
r in each positive direction (e.g., ft = x,y,z in 3D). It 
is straightforward to show that equations of motion for 
phase variable Or are giV|^n by Eq. (|6|) with the substitu- 
tion H given by Eq. ([? 



dOr 

'dt 



SH[Or, A] 
SOr'it) 



+ Cr. 



(8) 



In order to get a closed set of equations we further have 
to specify the dynamics of the twist variables A^, which 
is simply the average phase difference between opposite 
faces on the d-dimensional hypercube. In the absence of 
external currents, the physical boundary condition, cor- 
responding to an open boundary in real systems, should 
satisfy the condition that-there be no current across the 
boundary, which leads tola 



or, equivalently, 
dA„ 



dt 



SH[Or,A] . 



(9) 



(10) 



where Fa = l/L"^. As shown in Ref. || the noise term 
satisfies iC^it)) = {C^{t)Cr{t')) ^ and {C^{t)Cm = 
2TTASfj,uS(t). We term the dynamics defined in this way 
[by Eqs. (|) and (0)] RSJD with the FTBC, whereas 
the RSJD with the PBC is given by Eq. (|) with H in 
Eq. d). 



2. Relaxational dynamics 

Next we introduce the simpler phenomcnological re- 
laxational dynamics called time-dependent Ginzburg- 
Landau-Langevin-dynamics, which represents a noncon- 
served dynamics,E2l for the complex order parameter ipr 
on a discrete lattice: 



d'ipr 



SMt) 



(11) 



where F is the diffusion constant, Fgl is the discrete 
version of the GL free energy functional (H), and the 
white noise term satisfies (Cr(0) = and (Cr(*)Cr' (0)) — 
2kBTr6rr'S{t). The order parameter relaxes towards a 
configuration which locally minimizes the free energy, 
and the noises force the metastable states to decay. In 
the London limit the system can be described solely by 
the phase Or{t) of the order parameter ipr = '0oe'^'' with 
■00 fixed to a constant. Hence, by neglecting the am- 
plitude fluctuations and discretizing the time-dependent 
Ginzburg-Landau equation of motion, we find the phase 
equations of motion for the RD defined by 



dOr 

1t 



SH\ 



SOrit) 



(12) 



where H is the XY Hamiltonian (|l]) in units of J, the 
time unit is h/TJ, and the dimensionless thermal noises 
satisfy (Cr(i)) = and 



(Cr(t)Cr'(O)) =2TJ(t)5r 



(13) 



with T in units of J/ks- From Eq. ([121), the RD equa- 
tions for the phases in the case of the PBC are given 

by 



dOr 
1t 



E 



sm 



fi'r+fi) + Cr, 



(14) 



with periodicity on the phase variables: Or — Or+Lp.- 

We now proceed to the case of the FTBC for RD. 
In this case, in addition to the equations of motion for 
phases (Eq. (|2|) with substitution H[Or\ by ^^[^r. A] in 
Eq. (|^) ) we need dynamic equations for the twist vari- 
ables A^. Relaxational dynamics means that these equa- 
tions are of the form 



dA^ 

dt 



= -F 



5 H[0r,A] 

<5A, 



A- 



+ c 



(15) 
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which is identical to the form derived for RSJD [see 
Eq. ( |l0| ) , where Fa — 1/ L*^ was determined from the re- 
quirement that no current flows through the boundary] . 
We here define the dynamic equations for in the RD 
case with the same value of Fa , which makes the equa- 
tions identical to the corresponding equations in RSJD. 
Within the same interpretation that /rr+/i ~ 4>rr+ii and 

-^rr-l-A ^ Sin((/)rr+£)_with 0rr+A = Or - ^'r' " aS for 

RSJD (see Sec. [I C 1 ), we are again imposing a condi- 
tion consistent with that there be no current across the 
boundary. 

In the simulations, the coupled equations of motion 
are discretized in time (we use the discrete time step 
At = 0.05 and 0.01 for RSJD and RD, respectively) and 
numerically integrated using the second order Runge- 
Kutta-Helfand-Greenside (RKHG) algorithm,E£l which is 
much more efficient than the first-order Euler algorithm 
since it can reduce the effective temperature shiftO due to 
the discrete time step significantly. In the case of RSJD 
we apply the efficient fast Fourier transformation method 
(see, for example, Ref. 0), which makes the overall com- 
puting time 0{L'^log2 L) in d dimensions. [For compar- 
ison, the RD simulation requires 0{L'^).] The thermal 
noises are generated from a uniform distribution, whose 
width is determined to satisfy the noise correlation (see 
above) at a given temperature. 



III. SCALING RELATIONS 

A. Scaling in equilibrium 

In order to obtain the dynamic critical exponent z from 
equilibrium fluctuations of the system we use two differ- 
ent scaling relations: One is the flnite-size scaling of the 
time correlations of the supercurrent and the other is t 
finite-size scaling of the linear resistance. Fisher et ai, 
proposed a general scaling theory of the conductivity for 
a homogeneous superconductor, which has_been studied 
further explicitly by Dorsey and co-worker oll3 The pre- 
dictions from this scaling theory are very general and 
depend only on the dynamic scaling assumption and the 
existence of a diverging correlation length ^ |T — Td 
From a simple dimensional analysis, it is easy to show 
that the order parameter scales as ^ ^ ^j^^g 
the superfluid density scales as ps ~ [V'l^ Below 
Tc one has a^tu) ~ ips/i^, and accordingly one deduces 
thai, the frequency-dependent linear conductivity scales 
aS 



(16) 



where is a universal scaling function, the dynamic 
critical exponent z is introduced from t , and r is 
the characteristic time scale. Precisely at Tc, Eq. ( p^ ) 
turns into the finite-size scaling form of the conductivity: 



This scaling relation can be put to practical use in the 
case of the PBC because for this boundary condition ps 
has the required size scaling. On the other hand, it can- 
not be used for an open boundary condition or foe, the 
FTBC because in these cases ps = at any L and T.El For 
the FTBC we will then instead use the finite-size scaling 
of the linear resistance described below. 



1. Scaling of supercurrent correlations 

The conductivity cr(w) may be related to the supercur- 
rent correlation function G{t), which for the XY model 
in d dimensions is given by 



G{t) = {F{t)Fm , 



(18) 



where the global supercurrent F{t) flowing in a given di- 
rection, say, X, is written as 



^^(i) =^sin(0r-er+x). 



(19) 



The correlation function G{t) is a key quantity in de- 
scribing the dynamic response of vortex fluctuations and 
is for t = directly related to the static helicity modu- 
lus.Q The connectioa. between cr{uj) and G{t) in the RSJD 
case is expressed asQ 



a{uj) = 1 + 



iPs 



1 

TJo 



dte"^'G{t), 



(20) 



where the conductivity is measured in units such that the 
shunt resistance _Ro = Ij and the superfluid density ps is 
given by 



Po(l-^G(t^O) 



(21) 



with the bare superfluid density po = (cos(0r — (^r+x))- 
The dynamic dielectric function l/e(a;) in 2D is also ex- 
pressed aal3 



Re 



Im 



1 




1 LU 








1 




[e(c^)J 


PoT Jo 



dtsmLotG(t), (22) 



(23) 



where 



(24) 



(17) 



The helicity modulus 7 corresponds to the superfluid den- 
sity Ps and is given by 7 = ps = po/e{0). The conduc- 
tivitv, cr(w) in RSJD can be further siniplifled into the 
formQ 
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cr(w) = 1 - 



1 Po 

e{uj) 



(25) 



Expressing the scaling in terms of G{t) leads to the scal- 
ing form 



G{t) = e-'^FGitrn. 



(26) 



which at Tc for 3D turns into the finite-size scaling form 
(see Appendix ^) 



LG{t) = FcitL-^), 



(27) 



while in 2D a logarithmic correction (see Appendix ^) 
needs to be included 



In ( - ) G{t) = FaitL- 



(28) 



where Fq{x) is the scaling function for G{t). In the fol- 
lowing we will use the scaling relations Eqs. (|2^) and (|27| 
in 3D with the PBC and Eq. (pi) in 2D with the PBC. 



2. Resistance scaling 

In order to obtain a finite-size scaling at criticality for 
the FTBC for which, like for any open boundary condi- 
tion, ps — at any temperature and any lattice size, we 
relate the resistance R to the fluctuations of the twist 
over the sample. The voltage across the sample in the /t 



direction Vn 



-LAfj_ (see Ref. ph and the linear resis- 



tance Rfj^ in the same direction are related to the voltage 
fluctuation by the fluctuation-dissipation theoremll3 



Rn 



2T 
1 



2TQ 



dt{v,,{t)v^m 



([A^(e)-A^(0)]2) 



(29) 
(30) 



where the approximation becomes exact for a sufficiently 
large time Q, as shown in Appendix ^ (a similar approx- 
imation has been used for RSJD with open boundary 
condition in Ref. p^ . In the present simulation we use 
O = 2000 and perform average over all d directions, i.e.. 

Since i?^ scales as the inverse of the characteristic time 
scale-jji the critical region, the finite-size scaling takes the 
formca 



R=1-Fr{{T~T,)L'I''), 



(31) 



where v is the critical exponent for correlation length 
^ |r — Tel") and Fji{x) is the scaling function for 
R. Precisely at T^, Ffi{x) = -F/?,(0) becomes a constant 
independent of L and we get 



R^ L' 



(32) 



which can be used to determine z, once Tc is known. The 
resistance scaling can also be turned into an intersection 
metho(S3 for determining z and using that 



\n{RL/RL' 
ln{L/L') 



In 



Fn{{T - Tc)L^/-)/Fn{{T - Tc)L'^''') 



\tl{L/L') 



(33) 



for two different lattice sizes L,L' . Thus, if we plot 
\n{RL / Rl') / \n{L / L') as a function of temperature for 
several pairs of sizes {L-L'), all curves intersect at a single 
unique point (Tc,— z).E3 Once Tc and z are determined 
through the above intersection method, all data can be 
made to collapse onto a single scaling curve by plotting 
RL^ as a function of the scaling variable (T — Tc)L^^'^ 
with the correct value of the exponent v [see Eq. (pi|)]. 



Scaling of relaxation towards equilibrium: 
Short-time relaxation 



Recently, it has been found that a universal scaling in 
time can also be constructed for the relaxation towards 
equilibrium when starting from a nonequilibrium config- 
uration. Since such a relaxation is usually rather fast,rit| 
is often referred as the short-time relaxation method.EHl 



By this method several critical exponents have been suc- 
cessfully determined for the unfrusteated and the fully 
frustrated Josfphson junction arraynJ as well as for the 
Ising modeLcjS In these studies Glauber dynamics in 
MC simulations has been used to obtain time series of 
measured quantities, such as the magnetization and the 
Binder's cumulant. Here we apply this method to the 
XY models with more re alisti c dynamics, RSJD and RD, 
both introduced in Sec. [IC, in order to determine the 
value of thejdvnamic critical exponent z. For convenience 
we measurec3 



sign 



cos 9y-{t) 



(34) 



starting from the initial condition 0r(O) = 0. Since 
'0(i = 0) = 1 at any system size L, the finite-size scaling 
form becomes 
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^{t,T,L) = F^{t/LMT-T,)L'/''), 



(35) 



with the scahng function Ftp{x, y) depending on two scal- 
ing variables, satisfying (0, y) = 1 at any y. At Tc, z is 
easily determined from Eq. ( pq ) because in this case the 
second argument of the scaling function vanishes and ip{t) 
curves obtained for different sizes can be collapsed onto a 
single curve when plotted against the variable t/L^ . We 
can also determine Tc by an intersection method similar 
to Eq. ( p3| ) as follows: If the first argument of the scal- 
ing function is fixed to a constant {t/ L'^ = a) for a given 
system size L and z, then has only one scaling variable 
(T — Tc)L^^'^ , and can thus be written as 



(36) 



Accordingly, if we plot ^ with fixed a as a function of T 
for various L, all curves should intersect at Tc. However, 
because a depends on the value of z which cannot be 
independently determined by this method we start the 
intersection method from the z value determined from 
the scaling collapse at Tc- The values of Tc and z ob- 
tained from this intersection method can be refined by 
the iterating intersection construction. Finally, to exam- 
ine the consistency we collapse the data for all tempera- 
tures and lattice sizes onto a single scaling curve in the 
variable {T — Tc)L^^'^ at fixed a = t/L^ , which in addition 
is a check of the consistency against the known value of 
the static exponent v. 



IV. SIMULATION RESULTS 

A. 2D XY model 

In two dimensions, there has been some contro- 
versy over the value of the dynamic critical expo- 
nent: There has been a theoretical aaproach by Ambe- 
gaokar, Halperin, Nelson, and SiggiacJ (AHNS) predict- 
ing ZAHNS = l/2er^'^, where the Coulomb gas (CG) 
tem peratur e T^^; = T/2ttpo and 1/e = l/e(0) (see 
Sec. [II A 1 ). On the other hand, a si»ple scaling argu- 
ment has yielded Zscaic = l/eT*"^*^ — 2.0 Also, in numeri- 
cal simulations, there have been some differences: On the 



one hand, zahns has been observed in Ref. 25 from RSJ 
simulations, while ZscaisL has been concluded for RSJD 
and RD (Refs. ^ and 2^, for Langevin dynamics of CG 
gas particles (Ref. E6), and for the MC simulation of 
lattice CG (Ref. |l9|). Although the question is not com- 
pletely resolved yet, we strongly believe that when the 
fluctuating twist boundary conditionO (see Ref. ^ for a 
comparison between a conventional boundary condition 
and the FTBC) is used, ^scaic is the correct result. Al- 
though the above mentioned two z values are different 
below the KT transition, they give the same value of 2 
at the KT transition. In Ref. 28, however, z « 1 was 



PBC, while in Refs. |2| and |^ a very large value z « 5 
has been suggested from a scaling analysis of existing ex- 
perimental data and from an analytic calculation using 
Mori's technique, respectively. 

In the low-temperature phase of the 2D XY model, we 
can alternatively derive Zscaic in the following way: The 
potential barrier, which a bound vortex-antivortex pair 
should overcome in order to escape, is given by 



AV 



T 



InL, 



and the escape ratio T ^ exp(— AV^/T) for one pair is 
simply related to the total probability of escape, P, by 

P = L^nV, 

where n is the vortex pair density. The time scale r of the 
system is inversely proportional to P and thus is given 

by 



exp(Ay/T) 



and we obtain the dynamic critical exponent 



1 



gyCG 



concluded from a simulation of RSJ dynamics with the 



in accordance with Ref. where z has been obtained 
from a simple scaling argument and the observed 1/t be- 
havior of the correlation function G{t). 

In this section, we investigate the dynamic critical ex- 
ponent of the 2D XY model with RSJD and RD at, be- 
low, and above the KT transition. We use the FTBC as 
well as the conventional PBC and use various methods 
such as the resistance scaling, the scaling of the supercur- 
rent correlation function, and the short-time relaxation 
method. The results are summarized in Table |[ As 
seen from Table || only the FTBC gives results in accor- 
dance with the [expected value: Zscaic = -zahns ~ 2 at 
T = 0.90(w T,),y whereas z^caic « 3.4 and zahns « 2.8 
at T = 0.80.B Furthermore, this is the case both for RD 
and RSJD. In contrast, the results for the PBC are in- 
consistent both with Zscaie and zahns ■ From this we con- 
clude that the FTBC is an adequate boundary condition 
in the context of open systems like superfluid and su- 
perconducting films. It is also interesting to note that 
also the short-time relaxation method for RSJD with the 
FTBC gives results consistent with Zscaic- The results in 
Table || will be further discussed in Sec. 0. In the follow- 
ing we present the simulation results on which Tabic | is 
based. 



1. Critical temperature 

First we fix the temperature to T = 0.90 « Tc and 
focus on the dynamic critical behaviors at the KT tran- 
sition. The results from the resistance scaling R oc 



6 



for the FTBC (see Sec. |III A 2| ) are displayed in Fig. | (the 
data points are taken from Ref. |8|), where the slopes of 
the lines in the log-log plot correspond to z w 2 for both 
RSJD and RD. Consequently, our result z w 2.0 is i n ac- . 
cordance with other existing theoretical predictiong23'Efl 
while it contradicts recently suggested, very large values 
in Refs. |29| and 

In order to determine z at the KT transition for the 
PBC we use the following finite-size scaling form Eq. ( |2^ ) 
(see Sec. [II A 1) of the supercurrent correlation function 
G{t) 

In G{t) - FaitL-'). 

Figure ^ shows the corresponding scaling plot at T = 
0.90 w Tc both for RSJD and RD [Figs. |(a) and (b), 
respectively]. Very good scaling collapses are obtained 
in both cases with z — 1.5 for RSJD and z — 2.0 for 
RD. This clearly demonstrates that the value of z for 
RSJD with the PBC is different from the expected value 
of 2 which was obtained with the FTBC. In these scal- 
ing collapses one should note that the relaxation is much 
faster for RD than for RSJD, as is apparent by compar- 
ing the scales on the horizontal axes [note that vertical 
axis is in a logarithmic scale in Fig. ^(a) and in a linear 
scale in Fig. 0(b)]. It is also interesting to note that the 
value of the constant c in the logarithmic correction of 
Eq. ( p8| ) comes from the static properties, as described 
in Appendix and consequently should be independent 
of the dynamics. In accordance with this expectation the 
good scalings in Fig. ^ are achieved with the same value 
of c: Both for RSJD in Fig. |(a) and for RD in Fig. |(b), 

we found that c — 0.60 gives a good collapse. 

In Fig. H, we next show the decay of -0 (see Sec. [II B 
for details) at T = 0.90 « for RSJD with the(a) 
FTBC and (b) PBC, which demonstrates that z « 2.0 
(for the FTBC) and z « 1.2 (for the PBC) result in 
good data collapses to scaling curves. However, only the 
FTBC leads to the expected value z « 2.0. One should 
also note that the PBC results in an extremely slow de- 
cay of "i/j. Similarly Fig. ^ shows the decay of 'ip for RD at 
T = 0.90 with the (a) FTBC and (b) PBC. In both cases 
good data collapses are obtained for z w 2.0. In this case 
of RD, both boundary conditions have the same magni- 
tude of the decay time scale. A possible interpretation is 
discussed in Sec. 0. 



Figure || shows the finite-size scaling of resistance at 
T = 0.8(< Tc) for 2D RSJD and RD with the FTBC 
(all data are from Ref. H) and we find z « 3.3 for both 
types of dynamics. In Ref. ||, this value has been com- 
pared with Zscaie and zahns at this temperature and it 
has been concluded that the observed value 3.3 is very 

close to Zscalo « 3.4. 

In the same way as in the previous section the 
temperature-dependent z in the low-temperature phase 
can also be probed by the short-time relaxation method 
described in Sec. IIIB. The divergence of the correlation 
length in the whole low-temperature phase in 2D turns 
the finite-size scaling form (p5[) into the simpler form 



(37) 



i; = F^{t/L'), 



with the temperature-dependent z. Figures I and 
show the finite-size scaling of the short-time relaxation 
at T = 0.80 with the (a) FTBC and (b) PBC. The value 
z « 3.2 found in Fig. |(a) for RSJD with the FTBC is 
in agreement with z ~ 3.3 obtained from the resistance 
scaling in Fig. ^ within numerical accuracy. We inter- 
pret this as an evidence that z for the RSJD with the 
FTBC is indeed given by Zscaic- This is in contrast to 
the RSJD with the PBC for which at the same temper- 
ature (T = 0.80) z « 1.4 is determined from short-time 
relaxation as shown in Fig. ^(b). As will be discussed in 
Sec. we interpret this as further evidence that, in case 
of the 2D RSJD, z does depend on the boundary con- 
dition. The short-time relaxation for RD gives a quite 
different result: z « 2 is obtained for T = 0.80 for both 
the FTBC and PBC, as shown in Fig. @. If one compares 
this with the results at Tc in Fig. ^, where z « 2 is also 
obtained, the implication is that the result z « 2 for the 
short-time relaxation is expected at any temperature in 
the low-temperature phase both for the FTBC and PBC. 
As will be discussed further in Sec. this suggests that 
the short-time relaxation for RD does not probe the true 
equilibrium critical dynamics. 



3. High-temperature phase 

In the high-temperature phase there is a finite screen- 
ing length ^ which diverges as Tc is approached from 
above. Close to Tc one then expects that the character- 
istic time scales as 



2. Low-temperature phase 

The 2D XY model is special in that the whole low- 
temperature phase is "quasi" critical. This means that 
each temperature in the low-temperature phase is char- 
acterized by a temperature-dependent dynamic criti- 
cal exponent z. Just as in the previous section, this 
temperature-dependent z can be determined from the 
size scaling of the linear resistance, i.e., R oc L^^. 



In case of the PBC, we can estimate ^ and t following 
the method in Ref. ^ ^ is obtained from the wavevector 
dependence of the static dielectric function l/e(0) intro- 
duced in Eq. (|2^ ) . The characteristic frequency ujq ^ 1/t 
is determined from the frequency dependence of l/e(a;); 
luq is the position of the dissipation peak in |Iml/e(Li;)|. 
The resuh for RSJD with the PBC is shown in Fig. |, 
where z « 2 is found from o^o ^ C^^. It should be noted 



7 



that since this result is obtained for a temperature range 
where ^/L <C 1 it is expected to be independent of the 
treatment of the boundary and hence applies to both the 
PBC and FTBC. The same method apphed to RD also 
gives z « 2.0 as shown in Ref. In Fig. || the dotted 
line with slope —1 is also shown and corresponds to the 
result in Ref. |2^, where z « 1 was obtained for RSJD 
with the PBC in the same temperature range. Conse- 
quently, Fig. plimphes, in contrast to Ref. 28, that z = 1 



is the correct value for the PBC as well as for the FTBC, 
when z is determined from r ~ 



B. 3D XY model 

Next we turn to the 3D XY model with current con- 
serving RSJD and nonconserving RD, respectively. Both 
dynamic models have been used to describe the dynamic 
properties of high superconductors. Whereas it is gen- 
erally agreed upon that the static critical properties are 
those of the 3D XY model in a region close to with 
the corresponding static critical exponents,E3 there is less 
consensus on the dynamic critical propertiesj— ■ Several 
seemingly rnutualLy inconsistent experimentalE^ED and 
simulationalEjH'EJ results have been r eport ed. Similarly 



to the 2D case described above in Sec. IV A, we will here 



arrive at a somewhat entangled picture by comparing 
values of z obtained from the scalings in equilibrium for 
the two types of dynamics (RSJD and RD) with the two 
types of boundary conditions (the PBC and FTBC), as 
well as from the short-time relaxation method by observ- 
ing the time evolution towards equilibrium when starting 
from a nonequilibrium configuration. For convenience, 
the results of the simulations for the 3D XY model are 
summarized in Table O. 



1. Resistance scaling 

We start with the determination of z for the FTBC us- 
ing the finite-size scaling of the linear resistance, which 
is calculated from the equilibrium fluctuation of twist 
variable A [see Eq. (|30|)] A shorter presentation of these 
results has also been given in Ref. In 3D the cor- 
relation length diverges as ^ ^ \T — T^^"^, making the 
extended scaling form of Eq. (pl|), as well as the inter- 
section method in Eq. (33) applicable in addition to the 
relation R ^ at Tc- 

We first present the result for the scaling of the linear 
resistance for RSJD with the FTB C. By using the in- 
tersection method explained in Sec. [II A 2 [see Eq. (33)] 
we determine Tc and z simultaneously from the unique 
intersection point, as shown in the inset of Fig. ^, which 
yields ^ 2.200 and z « 1.46. We then display in the 
main part of Fig. ^jthe scaling plot of the linear resistance 
[see Eq. (pl|)], RL^ as a function of the scaling variable 
r7> 



(T - Tc)l7''', at T = 2.17, 2.19, 2.20, 2.21, and 2.23 for 



L = 4, 8, and 16. Here we used z and Tc found from 
the intersection method, i.e., z — 1.46 and Tc = 2.200, 
respectivelv, and the known value of the static exponent 
v ~ 0.67 S3 We also tried to vary the values of Tc, z, and 
ly in the scaling plot and concluded that z = 1.46 ± 0.06 
for the case of 3D RSJD with the FTBC. It is notewor- 
thy that Tc « 2.200 from the intersection method is very 
close to the known value of Tc [w 2.2018 (Ref. |l|)] from 
the MC simulation. 

For RD with the FTBC we only focus on the scaling re- 
lation R ^ L^^ at Tc, since it is found that the resistance 
for the RD case is harder to converge due to a sensitiv- 
ity of the result on the discrete time step in the numer- 
ical integration of dynamic equaticms even when using 
the second-order RKHG algorithm.Ej In contrast, we did 
not observe any significant sensitivity to the time step in 
RSJD and we fix At — 0.05 throughout the present work 
for RSJD. In order to overcome the problem in RD due 
to the finite-time step we obtain data for two different 
time steps {At = 0.05, 0.01) and linearly extrapolate to 
At = 0, as shown in Fig. |(b) for L = 4, 6, 8, 10, 12, and 
16. The slope of the line in the log-log plot of R versus L 
in Fig. |(b) gives 2 « 1.5 also for 3D RD with the FTBC. 



2. Supercurrent scaling 

For the PBC we use the scaling of the supercurrent cor- 
relation function G{t) introduced in Sec. IIIAl in order 
to obtain z. In Fig. |l^(a) we use the finite-size scahng 
form in Eq. ( |27| ) and plot LG as a function of the scal- 
ing variable t/L"" for (a) L = 8, 16, and 24 for RSJD 
and (b) L = 6, 8, 12, 16, and 24 for RD, respectively: 
Optimal data collapse is achieved for z = 1.5 (RSJD), 
and z — 2.05 (RD), respectively. In the insets of Fig. |l^ 
we use (a) z = 2.0 for RSJD and (b) z = 1.5 for RD, 
respectively, and show that the data collapse becomes 
significantly worse and consequently conclude that the 
z values obtained by this data collapse method are well 
determined [see the main parts of Figs. [0|(a) and|l^(b)]. 
One notes that for RSJD z w 1.5 is obtained for both the 
FTBC and PBC, whereas for RD z « 1.5 and z « 2 are 
obtained for the FTBC and PBC, respectively. 

In the critical region above Tc where ^ <C L we in- 
stead have £,G{t) as a scaling function with the scaling 
variable t/£,^ [see Eq. (p6|)]. Figure |ll| shows this scal- 
ing results at temperatures above Tc (T = 2.25, 2.30, 
and 2.40) for L = 24. By comparing with the results 
for L = 32, it is explicitly checked that there remain 
no significant finite-size effects in the current tempera- 
ture range. In Fig. ^ the corresponding iialues of ^ are 
taken from high precision MC simulations Jld As shown in 
Fig. the optimal value z = 1.4(2) is found for RSJD 
and z = 1-9(2) for RD, respectively, which is consistent 
with the finite-size scaling of G{t) at Tc. However, we 
note that the determination of z in the case of the finite- 
^ scaling above Tc yields a somewhat larger uncertainty. 
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Furthermore, z « 2 found for RD with the PBC is par- 
ticularly intriguing to understand since we expect that 
this result should be independent of boundary condition 
in this high-temperature regime where ^ <C i: Thus one 
expects the same value z « 2 for the FTBC in this high- 
temperature regime. This in turn suggests a discontinu- 
ous jump in the z value from z « 2 to z f» 1.5 as Tc is 
approached from above, since z w 1.5 at Tc was observed 
in the scaling of the linear resistance in Sec. IVB 1. This 
possibility is also discussed in Sec. 0. 



3. Short-time relaxation scaling 



The short-time relaxation method described in 



Sec. [II B probes the relaxation towards equilibrium from 
a nonequilibrium configuration. We start with the pre- 
sentation of the results obtained for RSJD and RD with 
the FTBC. Using the scaling form of -0 in Eq. (H) at Tc, 
where the scaling function has only one scaling variable 
t/i^, we first show in Fig. ^ the scaling plot of ■)/' at 
T = 2.20 for (a) RSJD with L = 4, 8, and 16 and (b) 
RD with T = 6, 8, and 10, respectively. All the data 
can be made to collapse onto a single curve in a broad 
range of the scaling variable for z « 1.5 and z « 2.0 for 
RSJD and RD, respectively. However, the above method 
presumes a priori knowledge of Tc. To circumvent this, 
one can alternatively use an intersection method with 
a fixed value of a = tL~^ in the fi rst a rgument of the 
scaling form in Eq. ( ^6| ) (see Sec. IIIB). In insets of 
Fig. P we display data points at T = 2.17, 2.19, 2.20, 
2.21, and 2.23 for (a) RSJD with L 4, 6, 8, and 16 and 
(b) RD with T = 4, 6, 8, and 10, and show the results 
of the iterative intersection method. We obtain again 
z ~ 1.5 and z « 2.0, as well as the estimations of the 
critical temperatures Tc « 2.200 and Tc « 2.194 for (a) 
RSJD and (b) RD, respectively. We believe that the ex- 
istence of an unique intersection point in each dynamics 
with the value Tc ~ 2.200 obtained for RSJD, which is 
in very good agreement with Tc « 2.200 obtained previ- 
ously from the resistapfie scaling, and with Tc « 2.2018 
from MC simulations,El make this short-time relaxation 
method very reliable. One notes that the slight temper- 
ature shift for RD is again the effect of the finite time 
step, as already observed in the calculation of the linear 
resistance. We have also checked the dependence on a 
values and observed no significant changes in resulting 
values of Tc and z in a broad range where 0.4 < ip < 0.9. 
Using z and Tc found from the intersection method, we 
in Fig. |l3| confirm that the full scaling form is borne out 
to high prceision with u = 0.67 determined from MC 
simulations .Eil 

We next consider the short-time relaxation for PBC, 
and show in Fig. n3 the scaling plot at T 2.20, 
ijj = F.^{t/L''), for (a) RSJD and (b) RD. Treating z as 
a free parameter, we obtain z ^ 1.2 and z « 2 for RSJD 
and RD, respectively. This suggests that the value for 



RSJD with the PBC is lower than z « 1.5 obtained from 
the same short-time relaxation method for RSJD with 
the FTBC, whereas for RD a value z « 2 is obtained 
both for the FTBC and PBC. As already observed in 
2D, RSJD with the PBC has a very large decay time 
scale in contrast to RSJD with the FTBC as well as to 
RD with both the PBC and FTBC. 



C. 4D model 

For completeness we also determine z in 4D.Ei As 
a prerequisite we first estimate Tc through the use of 
MC simulations in conjunction with the finite-size|-Sjcal- 
ing analysis of the Binder's fourth-order cumulantcJ U, 
which is independent of L precisely at Tc, 



U{L,T)^1 



3(|to|2) 



2\2 ' 



with the order parameter m = 'Y^^.e^^" / L"^ . The results 
are shown in Fig. |l^ and Tc w 3.31 is found from MC 
simulation pidth the PBC, which is consistent with ear- 
lier reportso but has a higher accuracy. From the MC 
simulations we also verified that v in 4D has the expected 
mean- field value v = 1/2 (see, e.g., Ref. ^). Since, as 
noted in the previous section, the size of the discrete time 
step in the integration of the dynamic equations of mo- 
tion can lead to an effective increase of temperature, we 
explicitly determine the effective Tc for RSJD and RD 
with the time step At = 0.05 from the crossing point of 
U{L,T): Figure |l^ shows that there is no significant dif- 
ference between the effective and nominal temperature 
for RSJD, leading to Tc(RSJD) « Tc(MC) « 3.31. On 
the other hand, for RD we from the crossing point obtain 
Tc « 3.25 at the same time step At — 0.05, in parallel 
with what was found for RD in 3D. It is to be noted that 
although the above critical temperatures have been ob- 
tained only with the PBC, the same critical temperature 
is expected also for the FTBC since all static quantities 
such as Tc should not depend on boundary conditions 
used. 

Once Tc is known from the calculation of the Binder's 
cumulant, we can use the simple finite-size scaling form 
Eq. ( |3^ ) for the linear resistance calculated with the 
FTBC by E g. M ) (we use 9 = 2000 for both RSJD and 
RD). In Fig. |16|(a), we plot the linear resistance R versus 
L at T = 3.31(RSJD) and T = 3.25 (RD), and from the 
least-square fit we find z « 2.1 for RSJD and z w 2 for 
RD, respectively. In addition, we also measure the short- 
time relaxation with the FTBC and present the result for 
RD at T = 6 and 8 in the inset of Fig. |6|(a) by using the 
simple scaling form at T = Tc = 3.25, i.e., = Fjp{tL~^)^ 
which yields z ss 2.0 in accordance with the result from 
the resistance scaling. For RSJD with the FTBC, we 
construct the intersection plot for the short-time relax- 
ation (similar to Fig. |l3| for 3D) as displayed in the in- 
set of Fig. |l|(b), and get z w 2.0 and Tc « 3.31 from 
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the unique crossing point. It is interesting to note that 
the critical temperature obtained here for RS JD with the 
FTBC is in a perfect agreement with Tc found from the 
Binder's cumulant method for the other boundary condi- 
tion, the PEC. We then make the fuU scahng plot for 'ip 
in the main part of Fig. ^(b) with the mean-field value 
V = 0.5 and the estimated values Tc ~ 3.31 and z — 2.0 
above, resulting in a very smooth collapse. 

In short, we get z 2 in 4D with the FTBC regardless 
of the dynamics we use (see Tabic ID for a summary of 
results); this is reassuring since the value z = 2 is usu- 
ally expected in 4D|-wliere the phase transition acquires 
a mean- field nature llil'EJ 



V. DISCUSSIONS AND COMPARISONS 

As is clear from the simulation results presented in 
Sec. 1^ for two-, three-, and four-dimensional XY mod- 
els, a very entangled picture emerges as regards to the 
dynamic critical exponent z. In this section we discuss 
the main features. 



A. Discussions of the 2D XY model 

We start our discussion with 2D (see Table | for sum- 
mary of results) and first focus on RSJD at the KT tran- 
sition. For a 2D supcrfluid and superconductor the most 
widely expected value is z = 2 although there have been 
a few different suggestions (Refs. p8|-^. The value z — 2 
can be inferred from the observed nonlinear current- 
voltage {I-V) exponent a = 3 (Refs. |,|4|j2^, and 
together with the scaling argument that a — z + iJli 
It may also be directly obtained from the simple ar- 
gument of the escape over the boundary presented in 
Sec. |IV A| with the result z = l/eT^^ - 2, combined 
I jump condition at the KT transition 
I which leads to z = 2 at the transition. 



with the unii 
1/eTC^G ^ 4J2 

For the 2D XY model, the KT transition temperature is 
Tc « 0.9 (Ref. |l|) and as seen from Table |, RSJD with 
the FTBC does give the expected value. However, RSJD 
with the conventional PBC does in fact not give the ex- 
pected value: the supercurrent scaling gives z w 1.5 and 
the short-time relaxation gives z w 1.2. 

In order to understand the role played by the bound- 
ary conditions we consider a system with an open bound- 
ary, which is appropriate to describe a superconducting 
film and a film of ^He in usual experiments. In such a 
case, when a vortex-antivortex pair is introduced into the 
ground state and then annihilated across the boundary, 
the system relaxes back to the original ground state. The 
FTBC has been designed to keep the advantage of the 
PBC, which reduces the finite-size effect compared to the 
free boundary condition, as much as possible, while al- 
lowing this relaxing back. This relaxijig back is, however, 
prohibited by the conventional PBCfl One may note that 



the escape-over-barrier argument in Sec. IV A implicitly 
presumes this relaxation back as a part of the escape 
process. One should also note that, when comparing to 
experiments with open boundaries, the FTBC has to be 
used in simulations instead of the PBC whenever the re- 
laxation process across the boundary is important. This 
perspective suggests that the observed difference between 
the FTBC and PBC at the KT transition for RSJD is due 
to the additional constraint on the physics caused by the 
PBC. 

This can be substantiated somewhat further by study- 
ing the low-temperature phase in 2D, where an ubiq- 
uitous "quasi" criticality with a diverging correlation 
length makes the critical finite-size scaling method ap- 
plicable. In Ref. H, z w 3.3 at T = 0.80 was found for the 
FTBC from the resistance scaling in agreement with the 
expected value z = 1/eT'-^'^ — 2 « 3.4 within munerical 
errors. However, an estimate of the equilibrium scaling 
for the PBC at T = 0.85 in Ref. |gave z « 1.6 instead of 
the expected result z ~ 2.8 (see Fig. 3 in Ref. ||). Thus in 
this low-temperature phase z determined with the PBC 
appears to be smaller (z < 2) than the one with the 
FTBC (z > 2). However, the value for the FTBC is the 
relevant one when comparing with experiments. 

The situation above Tc is as follows: The finite lin- 
ear resistance R calculated from the fluctuations of A 
for the FTBC [see Eq. (|30|)] can be related to the con- 
ductivity calculated for the PBC through the connec- 
tion R = Re[l/cr(a; = 0)] with cr(tj) in Eq. (|^. We 
have explicitly checked this relation in our simulations at 
T = 1.4, by comparing the two values for the FTBC and 
PBC, respectively, and found good agreement .Ha From 
this observation, we expect that in this high-temperature 
phase where the correlation length is smaller than the 
linear size of the system, R and z are independent of 
boundary conditions. Furthermore, in this temperature 
regime, transport properties like the linear resistance are 
dominated by free vortices (with density np) and accord- 
ingly we expect R oc np oc ^"^ (Ref. ||), leading to z = 2 
for both boundary conditions. However, from a computa- 
tional point of view, the calculation of the size-converged 
R for the FTBC in the high-temperature phase becomes 
difficult as we approach Tc from above, due to the diverg- 
ing correlation length. On the other hand, if we instead 
focus on the scaling of the characteristic frequency luq, 
which is expected to be proportional to R and can be 
calculated for the PBC, then we do indeed find an in- 
dication of the expected behavior, loq ~ seen in 
Fig. I for the PBC. For the FTBC this result is con- 
sistent with our observation z « 2.0 at Tc, whereas for 
the PBC the scaling at Tc gives z « 1.5, which differs 
from the expectation. Why is there then a difference in 
the PBC between z values at and above Tc? The point is 
that the long-time relaxation above Tc is governed by the 
thermally created free vortices, whose density satisfies 
np (X whereas the behavior precisely at Tc, where 
np = 0, is instead dominated by the bound pairs of vor- 
tices and antivortices. The conclusion is then that the 
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constraint imposed by the PBC on the vortex-antivortex 
escape gives rise to this pecuhar discontinuity of z pre- 
cisely at Tc. This is in contrast to the FTBC case where 
z appears to be a continuous function of T. 

Next we compare the results from the dynamic scal- 
ing in equilibrium and the short-time relaxation method 
which probes the relaxation when the system approaches 
equilibrium. For RSJD with the FTBC there is no differ- 
ence: the resistance scaling and the short-time relaxation 
method yield the same z at and below Tc (see Table |). 
However, for RSJD with the PBC the equilibrium scal- 
ing and the short-time relaxation scaling lead to different 
results, z w 1.5 and z « 1.2, respectively. In fact by com- 
paring Figs. |2|(a) and ^(b) one realizes that the approach 
to equilibrium from the chosen starting nonequilibrium 
configuration is much slower than the equilibrium relax- 
ation. Apparently the constraint imposed on the relax- 
ation by the PBC in combination with the nonequilib- 
rium starting configuration is causing the difference. 

We now turn to the discussions for RD, where for the 
FTBC we find from the resistance scaling the same z 
at and below as for RSJD (see Table Q). In this 
context it is interesting to note that the 2D XY model 
with the FTBC is dual to the lattice CG model with the 
PBC (see Ref. | for the mapping between two models), 
where the same values of the dynamic critical exponent 
{z = Zspajo = l/eT'^^ - 2) have been found in MC dy- 



namics.E3 Also, the continuum CG model with Langevin 
dynamics of the pure relaxational form has been found to 
give the same values of z.Ej Accordingly it is tempting to 
conclude that the result presented in this work for the 2D 
XY model with the FTBC is associated with the vortices 
and that it is essential to define the model so as to allow 
for a proper relaxation of vortex-antivortex annihilation 
across the boundary, which is not the case for the PBC. 
Furthermore, the result that z — Zgcaic appears to be uni- 
versal in the sense that it does not matter whether or not 
the underlying dynamics is purely relaxational (such as 
RD in this work, MC dynamics in Ref. |4^, and Langevin 
dynamics in Ref. p6| ), or it has an additional constraint 
like local current conservation in the RSJD case. 

Although the short-time relaxation method applied to 
RD gives the expected value z ss 2 at Tc, it fails to yield 
the equilibrium size scaling value below Tc- In addition, if 
we compare the decay behaviors at and below Tc, shown 
in Figs. ^ and respectively, we notice that the time 
scale of the relaxation does not depend significantly on 
the temperature or on the boundary condition, in sharp 
contrast to RSJD. We suggest the following reason: In 
RD the relaxations of spin waves and vortices are ef- 
fectively decoupled and the short-time relaxation in this 
case only probes the spin- wave degrees of freedom, which 
follow the purely relaxational dynamics with the trivial 
exponent z = 2 at any T, while the resistance scaling 
probes the vortex degrees of freedom. This is then in 
contrast to RSJD with the FTBC where both degrees 
of freedom are strongly coupled, leading to the same re- 
laxation time (and accordingly the same z value) for -ip 



and R. It is also interesting to note that z w 2 was also 
found in Ref. |l| from the MC simulations of the 2D XY 
model with the PBC at and below Tc by using a similar 
short-time relaxation method. 



B. Discussions of the 3D XY model 

We next turn to the 3D XY model (see Table |l| for a 
sum mary of results). The discussion for the 2D case in 
Sec. V A regarding the boundary conditions carries over 
to 3D, and we expect that the FTBC has to be used 
whenever the relaxation process associated with the ex- 
pansion and the subsequent annihilation of a vortex loop 
across the boundary is important because the conven- 
tional PBC prevents this relaxation. 

For RSJD, z w 1.5 is found from the linear resistance 
and the short-time relaxation method for the FTBC, as 
well as from the scaling of the supercurrent correlation 
at and above Tc for the PBC. In addition, the same value 
z « 1.5 is also found for RD with the FTBC from the 
finite-size scaling of the linear resistance. We note here 
that the MC simulations of the lattice vortex loop model 
in 3D (Ref. |l^) also have found the same value. The 
agreement between the z values for the three different 
dynamic models (RSJD and RD with the FTBC, and 
MC dynamics of the vortex loop model with the PBC) 
was also found in 2D. This value z « 1.5 obtained in 3D 
is consistent with z — d/2 (with d = 3 in 3D) for model 
E and model F describing critical dynamics of superfluid 
systems, iti-tJie classification scheme of Hohenberg and 
HalperinlB'Ej Consequently, it is again tempting to con- 
clude that the result for the 3D XY model can be associ- 
ated with the vortex loops and that the critical dynamics 
of RSJD and RD are equivalent as long as the boundary 
condition allows for the proper vortex loop escape over 
the boundary. 

As in 2D, we find that the short-time relaxation 
method for RSJD with the PBC (with the result z w 1.2) 
does not reflect the true equilibrium relaxation corre- 
sponding to z « 1.5, and we again suggest that this is 
due to the constraint imposed by the PBC. On the other 
hand, we find that the short-time relaxation method for 
RD with the FTBC gives z = 2.0(1) [see Figs. ||(b) and 
|l|(b)]. We propose the same explanation as we did for 
2D: The short-time relaxation in RD at criticality does 
not reflect the true long time relaxation because the vor- 
tex loop configurations are still out of equilibrium even 
when -0 w is reached. In this respect it is interesting 
to note that the values Tc = 2.20 and v = 0.67 used in 
the scaling collapse for ip in Fig. 13(b) with z — 2.0 agree 



with the value expected for the 3D XY model and that 
the same Tc — 2.20 was used in the resistance scaling in 
Fig. ||(b) and yielded z w 1.5. 

We next discuss the results for RD with the PBC. The 
scalings of the supercurrent correlation both at and above 
Tc (corresponding to the finite-L scaling and the finite-^ 
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scaling, respectively), as well as the short-time relaxation 
at Tc, consistently give z ~ 2. This value corresponds 
to model A of relaxational dynamics, in the Hohenberg- 
Halperin classification scheme. EJoEj The most striking 
feature in RD is that the scalings at Tc for the FTBC 
(resistance scaling) and PBC (supercurrent scaHng) cor- 
respond to different values, i.e., z ~ 1.5 and z « 2.0, 
respectively. 

In the high-temperature phase in 3D where ^ <^ L, o: 
expects that z is independent of boundary conditions. 
Consequently, z w 2 found for RD with the PBC at tem- 
peratures above implies z « 2 also for RD with the 
FTBC in the same hish-temperature regime, again con- 
sistent with model A.tll In contrast, z determined from 
the resistance scaling at Tc instead gives z « 1.5 for RD 
with the FTBC. We propose the same explanation for 
this discontinuity of z at Tc in the RD case with the 
FTBC as we did in 2D: Above Tc where f < i, the finite 
value of the resistivity reflects the physics of dissociated 
vortex loops whereas precisely at T^ where the resistiv- 
ity vanishes as L is increased, the physics is dominated 
by the large nondissociated vortex loops. 

We now compare our results in 3D with earlier studies. 
Values consistent with z « 1.5 have also been found in 
earlier simulations: z = 1.5(5) was obtained from the TV 
characteristics of the current-driven RSJ model with an 
open boundary (Ref. [52| ), and z — 1.5(1) was concluded 
from the scaling of the linear resistance for the MC sim- 
ulations of the XY model in the vortex representation 



with the PBC (Refs. |19| and 39) which corresponds to the 
FTBC in the phase representation as mentioned above 
and explained in Ref. |7|. Finally, MC spin dynamics ap- 
plied to the three component XY model gave z = 1.38(5) 
in Ref. On the other hand, the experinieptal situa- 
tion for high-Tc superconductors is less clear:E3 From sev- 
eral zero field dc conductivity experiments z « 1.5 has 
been found on single YBCO-123 crystalscfl and a similar 
result z — 1.6(1) was also obtained for a Bi-2212 crys- 
talHZl However, from the scaling of the magnetoconduc- 
tivity of a thick YBCO-123 film z — 1.25(5) was found in 
Ref. 36, whereas a similar experiment reported z w 2 in 
Ref. 55. From a theoretical point of view the renormal- 
ization group methods applied to the .relaxational model 
(model A) yield the resuh z = 2 + cr^Jlil'tj with 77 « 0.02 
and c « 0.7261, leading to z « 2.0. However, as far as 
we know, no corresponding calculation has been made 
for the 3D RSJ model. One may argue that since the 3D 
RSJ model is a bona fide model of a superconductor the 
critical dynamics should belong the dynamic uniijersal- 
ity class of model F which describes supcrfluidsO This 
gives z — d/2 = 1.5 for a model with the static properties 
given by the 3D XY modelB 



C. Discussions of the 4D XY model 

In case of the 4D XY model both resistance scaling at 
Tc and short-time relaxation give z « 2 for RSJD as well 
as for RD (see Table HI for summary of results) . This is 
in perfect accordance with the Hohenberg-Halperin clas- 
sification scheme where the RSJD case should be related 
to models E and F with z = d/2 = 2 and the RD case 
with the model A value z = 2. This in turn just reflects 
that 4D is the upper critical dimension. 



VI. SUMMARY 

We have found that the size scaling of the resistance for 
the XY model with the FTBC gives the dynamic criti- 
cal exponent z related to superfluid and superconducting 
systems with an open boundary. This is the case in two, 
three and four dimensions both for relaxational and RSJ 
dynamics. In 2D this applies for T < Tc, whereas in 3D 
and 4D the dynamics is critical only at T = Tc- However, 
the 3D case with relaxational dynamics has a discontinu- 
ity in the z value since the relaxation time r scales as 
T oc L^^^ at Tc, whereas it scales as t oc just above Tc- 

The short-time relaxation method, which probes the 
relaxation from a nonequilibrium configuration, does give 
the same result, except for the 3D case with relaxational 
dynamics where z w 2 is obtained. This discrepancy 
shows that although the short-time relaxation method 
very often is reliable and efficient, it cannot always be 
trusted jOS a determination of the critical equilibrium dy- 
namics.E^ 

The XY model with the PBC has a different dynam- 
ical size scaling behavior than with the FTBC. In 2D 
the XY model with the PBC and RSJ dynamics gives 
smaller values of z (z < 2) both at and below Tc- This 
demonstrates that the boundary condition influences the 
size scaling properties of the dynamics. Also in this case 
there is a discontinuity in z since t oc L^-^ at Tc but 
T cx just above. This is similar to the discontinuity 
found in 3D for relaxational dynamics. The short-time 
relaxation fails to give the equilibrium z at Tc for the 2D 
XY model with the PBC and RSJ dynamics. 

In 3D the XY model with RSJ dynamics and the PBC 
gives the same result as for the FTBC. Thus in this case 
the boundary condition does not influence the size scal- 
ing. This is in contrast to the 3D XY model with re- 
laxational dynamics which gives different result for the 
PBC and FTBC. 

In 4D all determinations give the simple relaxational 
value z « 2 independent of boundary condition and dy- 
namics. 

The actual values determined are consistent with the 
following sequences: the XY model with RSJ dynamics 
and the PBC at Tc is consistent with z = 1.5, 1.5, and 2 
for the 2D, 3D and 4D cases, whereas the FTBC is con- 
sistent with z = 2, 1.5, and 2. Similarly for relaxational 
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dynamics the PBC gives the sequence 2; = 2, 2, and 2 
whereas the FTBC gives z = 2, 1.5, and z — 2 for the 
2D, 3D, and 4D cases, respectively. 

The z values for a superconductor can be related to 
the nonlinear I^V.. f;xponcnt a through the scaling rela- 
tion a = 1 + z.Ejo The I-V measurements correspond 
to an open boundary and simulations with the FTBC 
are consistent with the scaling relation, as shown for the 
2D case in Ref. ||. On the other hand, the z values cal- 
culated with the PBC in 2D do not fulfill the relation 
because a > 1 + z for T < Tj,- In this sense, the bound- 



ary condition has direct physical significance. 
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APPENDIX A: SCALING FORM OF THE SUPERCURRENT CORRELATION FUNCTION 

For supercurrent correlation scaling in 2D, the superfluid density ps is proportional to the vortex dielectric function 
l/e(0) which in turn is plated to the conductivity by <j{Ld) oc l/[ia;e(w)]. Precisely at the KT transition l/e(0) has a 
logarithmic size scaling,EZ 

111 



eL(0) eoo(O) ln(L/c) 

, where c is a constant. This is consistent with the functional form in terms of a scaling function with a logarithmic 
correction for the frequency dependence 

11 1 , , 
77-r oc F^yujT). 



e(w) e(0) ln(i/c) 

The supercurrent correlation function G{t) is related to (7{lij) by a Fourier transform so that 



Im 



1 






-L 



dtuj cos ujtG{t) = 



\nL/c 



where F{x) — lm[F(a;)]. From this it follows that 

/CX3 1 poo 1 

^00 ^ J-00 

SO that ln{L/c)G{t) has the scaling form Fcit/T). 

In 3D we have instead 1/6^(0) oc ps oc 1/L at Tc and l/eoo(0) = so that this time l/eL(0) — l/eoo(0) oc 1/L. This 
means that going through the same steps as above for the 3D case gives the scaling form LG{t) = Fcit/T). 



APPENDIX B: APPROXIMATION MADE IN THE NYQUIST FORMULA FOR LINEAR RESISTANCE 



The step from Eq. ( |29| ) to Eq. (30) is equivalent to showing that 

dt{Amm ^ ^((A(e) - A(o))') 

in the limit of large Q. The left-hand side can, due to translational invariance, be written as 

r-e-s rO /.e 



1 /-^ . . If® 



s Je~s 



dt{A{t + s)A{s)). 



(Bl) 



The first term on the right-hand side is ((A(0) — A(0))^)/O, and the second term reduces to 
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e 



ds 



dt{A{t + s)A{s)) 



1 

e./o 



e 



ds([A(.)-A(0)]A(s)) 



(A2(e) - A2(0) - 2A(0)A(e) + 2A2(0)) 



26 



([A(e)-A(0)]^) 
29 



(B2) 

(B3) 
(B4) 



where we have used 2A(s)A(s) = dA /ds. Thus the sum of the first two terms on the right-hand side of Eq. (Bl) is 



equal to ([A(0) — A(0)] )/28 and it remains to prove that the third term vanishes in the limit 8 — > oo. This can be 
realized by changing the order of integration: 



dt(A{e + t + s)A{s)) = - 

_ 1 
~ e 

_ 1 
^ e 



dt / ds{A{Q + t)A{0)) 



dt{e + t){A{Q + t)A{0)) 



-e 
e 



dxx{A{x)A{0)). 



(B5) 
(B6) 
(B7) 



A finite relaxation time t means that (A(i)A(O)) oc exp(--t/T), which means that the last integral is finite for any 
finite T. This is the case at Tc whenever the system size is finite since r oc i^. Consequently, the third term vanished 
in the limit ^ oo for any finite L and the step from Eq. (^) to Eq. (|3^) follows. 
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TABLE I. Dynamic critical exponent 2 for 2D RSJD and 
RD with the PBC and FTBC. 
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at T = 0.90 










resistance scaling 
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supercurrent scaling 
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short-time relaxation 


1.2(1) 
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at T = 0.80 
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TABLE II. Dynamic critical exponent z for 3D RSJD and 
RD with the PBC and FTBC. 
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RD 




PBC 


FTBC 
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resistance scaling 




1.46(6) 




1.5(1) 


supercurrent scaling (L) 
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supercurrent scaling (^) 
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1.9(2) 




short-time relaxation 
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2.1(1) 


2.0(1) 


TABLE HI. Dynamic critical exponent z for 4D RSJD and 


RD both with the FTBC. 
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FIG. 1. Resistance scaling for 2D RSJD (open squares) and 
RD (solid squares) with the FTBC at T=0.9. The sohd lines 
represent R ~ L~'^ °. For both dynamics^ z ~ 2.0 is obtained. 
(The data points were taken from Ref. 
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FIG. 3. Short-time relaxation for 2D RSJD at T = 0.9 . 
(a) For the FTBC z « 2.0 and (b) for the PBC z ^ 1.2 are 
obtained. Note the enormous time scale for the PBC. 
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FIG. 2. Scaling of the supercurrent correlation function 
G{t) in 2D at the KT transition (T = 0.9) in case of the 
PBC for (a) RSJD and (b) RD. The scaling collapse of the 
data shown is for z = 1.5 and 2.0 in case of (a) RSJD and (b) 
RD, respectively. [The apparent spread of the data in (a) is 
only due to the logarithmic scale and insufficient convergence 
for the largest lattice sizes at large t.] 




FIG. 4. Short-time relaxation for 2D RD at T = 0.9 with 
the (a) FTBC and (b) PBC. For both boundary conditions 
2 ^ 2.0 is obtained from the data collapse. Note that the de- 
cays for both boundary conditions have the same time scale. 
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FIG. 5. Resistance scaling for 2D RSJD (open squares) and 
RD (solid squares) with the FTBC at T=0.8. For both dy- 
namics, z ~ 3.3 is obtained from the slopes of the lines in the 
figure. (The data points were taken from Ref. ^.) 



10 12 14 




FIG. 7. Short-time relaxation for 2D RD at T = 0.80 for 
the (a) FTBC and (b) PEC. For both boundary conditions, 
z ~ 2.0 is found from the best data collapse. 
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FIG. 6. Short-time relaxation for 2D RSJD at T = 0.80. 
(a) For the FTBC z « 3.2 and (b) for the PBC z x lA are 
obtained from the best data collapse. 
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FIG. 8. Characteristic frequency uio determined from the 
peak position of |Iml/e(a;)| vs the correlation length ^ for 2D 
RSJD with the PBC at temperatures T = 1.0, 1.1, 1.2, and 
1.3 (from right to left). The dynamic critical exponent z de- 
fined by cjo ~ is shown to have a value close to 2.0 (solid 
line) . For comparison, we also plot the dotted line which cor- 
responds to z — 1.0. 
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FIG. 9. (a) Scaling curve of the resistance R for 3D RSJD 
with the FTBC for L = 4, 8, and 16 at T = 2.17, 2.19, 2.20, 
2.21, and 2.23 with parameter values z — 1.46 and Tc = 2.200 
(both from the intersection method described in the text and 
shown in the inset), and the known value u = 0.67 (Ref. ^l| ). 
(b) Determination of z for 3D RD with the FTBC from the re- 
sistance scaling form R oc L^^ at Tc — 2.20. The data points 
are for L = 4, 8, 10, 12, and 16, and two integration time 
steps At = 0.05 and 0.01. Linear extrapolation to At = 
gives z — 1.47 from the least-squares fit. 



FIG. 10. Finite-size scaling at Tc of the correlation function 
LG{t) vs t/L^ for 3D (a) RSJD with the PEC and (b) RD 
with the PBC. In the main parts of (a) and (b), z — 1.5 and 
z — 2.05 are shown to give good scaling collapses for (a) RSJD 
and (b) RD, respectively, while in the insets the interchanged 
values [z = 2.0 and z — 1.5 for (a) and (b), respectively] are 
shown to be inconsistent with the scaling collapse. 
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FIG. 12. Short-time relaxations of ip in 3D with the FTBC 
a.t T = Tc — 2.20 are shown as functions of the scaUng vari- 
able tL"^ for (a) RSJD with L = 4, 6, and 16, and for (b) RD 
with L = 6, 8, and 10. From the scaling collapse z — 1.50 and 
z — 1.95 are found to yield smoothly collapsed single curves 
for (a) RSJD and (b) RD, respectively. 
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FIG. 11. Scaling of the correlation function ^G{t) above Tc 
is shown against the scaling variable for 3D (a) RSJD with 
the PBC and (b) RD with the PBC for L = 24 and T = 2.25, 
2.30, and 2.40. From the scaling collapse, with ^ from the MC 
simulation (Ref. ^l]), we obtain (a) z ~ 1.4 for RSJD and (b) 
z ^ 1.9 for RD. 
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FIG. 13. Finite-size scaling of the short-time relaxation of 
i> for 3D (a) RSJD for L = 4, 6, 8, and 16 and (b) RD for 
L = 4, 6, 8, and 10 both with the FTBC and at T = 2.17, 2.19, 
2.20, 2.21, and 2.23. As shown in insets, the scaling function 
i) = F^(t/L\ (T - Tc)L'^^'") with a fixed a = t/L^ suggests 
the existence of a single crossing point (a) at Tc = 2.200 
with {z,a) = (1.5,4.0) for RSJD and (b) at Tc = 2.194 with 
{z,a) = (2.0,0.5) for RD. All data points in the insets col- 
lapse onto a single smooth curve with the scaling variable 
(T -a^)L^/" for both models with the known value v = 0.67 
(Ref.Ej), as shown in main parts. 




FIG. 14. Short-time relaxation of il) in 3D with the PBC 
at T = Tc = 2.20 as functions of the scaling variable tL~^ 
for (a) RSJD with L = 4 and 8, and for (b) RD with L = 4, 
6, and 8. From the scaling collapse z « 1.2 and z ~ 2.1 are 
found for (a) RSJD and (b) RD, respectively. 
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FIG. 15. Binder's fourth order cumulant in 4D determined 
from MC simulations, RSJD, and RD for L = 4, 6, and 8 at 
several temperatures between T = 3.05 and T = 3.40. Data 
points for MC and RSJD coincide within error bars and give 
Tc ~ 3.31, whereas the RD results show a relative large dy- 
namic shift to Tc ~ 3.25. 
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FIG. 16. (a) Determination of z in 4D XY model with the 
FTBC from the resistance scahng form R ~ for L = 4, 
6, and 8 at Tc = 3.31 for RSJD and at Tc = 3.25 for RD (see 
Fig. |l^. From the slopes in log-log plot, z ~ 2 is concluded 
for both dynamics. Inset: Short-time relaxation of tp for 4D 
RD with the FTBC at T = Tc = 3.25 is shown against the 
scaling variable tL~^ for L — 6 and 8 and z = 2.0 is found, 
(b) Finite-size scaling of ii for 4D RSJD with the FTBC for 
L = 4, 6, and 8 at T = 3.25, 3.30, 3.31, 3.32, and 3.35. As 
shown in inset, the intersection method gives Tc ~ 3.31 with 
z = 2.0. The main part displays the full scaling plot of the 
form tl> = F^(t/L\ (T-Tc)L^''') with a = tL" = 2.8 and the 
mean-field value v = 1/2. 
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